A Lagrangian model for the evolution of turbulent magnetic and passive scalar fields 
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In this paper we present an extension of the Recent Fluid Deformation (RFD) closure introduced 
by Chevillard and Meneveau [l| which was developed for modeling the time evolution of Lagrangian 
fluctuations in incompressible Navier-Stokes turbulence. We apply the RFD closure to study the 
evolution of magnetic and passive scalar fluctuations. This comparison is especially interesting since 
the stretching term for the magnetic field and for the gradient of the passive scalar are similar but 
differ by a sign such that the effect of stretching and compression by the turbulent velocity field 
is reversed. Probability density functions (PDFs) of magnetic fluctuations and fluctuations of the 
gradient of the passive scalar obtained from the RFD closure are compared against PDFs obtained 
from direct numerical simulations. 

PACS numbers: 47.27.eb, 47.10.-g,02.50.Ey,52.30Cv, 52.35.Ra, 52.65.Kj 
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I. INTRODUCTION 



In recent years several Lagrangian type closures have 
been developed to model fluctuations of turbulent incom- 
pressible fields [J-Q. The starting point for all these 
models goes back to the so-called Restricted Euler clo- 
sure [a, [6| which models the time evolution of the gradi- 
ent tensor Aij := d%Uj 



dt 



I AikAki 



3 



k^kn 



(i) 



This model is derived using the incompressibility con- 
straint and neglects the viscous term and the anisotropic 
part of the pressure Hessian. The Restricted Euler model 
is already able to capture geometric features of vortex 
stretching and alignment but fails as a robust model to 
obtain stationary statistics due to the appearance of fi- 
nite time singularities. 

Several strategies have been applied to regularize the 
restricted Euler model where the most prominent ones 
are the Tetrad Model @ which models the anisotropic 
part of the pressure Hessian and the Recent Fluid De- 
formation approach which models both the viscous term 
and the anisotropic part of the pressure Hessian 

In this paper we extend the Recent Fluid Deformation 
approach to the case of the time evolution of the gradient 
of passive scalar and to kinematic magnetohydrodynamic 
turbulence. Comparing these two examples is especially 
interesting since both examples differ mainly by a sign in 
the stretching term. Thus structures which are expanded 
in one equation (e.g. gradient of passive scalar) are com- 
pressed in the other (e.g. MHD). Comparison with direct 
numerical simulations are performed to test the range of 
applicability of the extended Recent Fluid Deformation 
models. 



II. VELOCITY GRADIENTS 

Before we apply the Recent Fluid Deformation ap- 
proach to the gradient of a passive scalar and to kine- 
matic MHD turbulence, we first recall the main steps in 
the derivation of this model. Detailed information can be 
found in [l|, H|. Starting with the Navier-Stokes equation 

<9 f u + u- Vu + Vp = ^Au , V-u = (2) 

one takes the gradient 



d t d k Ui + d k (ujdjUi) + d ki p = vd kjj Ui 



(3) 



and denotes the velocity gradient diUj as Aij. Using the 
Lagrangian material derivative 4z = (dt +Ujdj) eqn. §5§ 
reduces to 



dt 



' Ai k A 



-dijp + vdkkAi 



(4) 



Assuming the pressure Hessian to be isotropic d^p = 

^f-dk k p and neglecting the viscous term would result 
in the Restricted Euler model given by eqn. ([1]). In 
order to improve the model for the pressure Hessian, 
Chevillard and Meneveau [l[ consider the pressure Hes- 

sian P nm = q X q X expressed in the Lagrangian frame 
and postulate that the Lagrangian pressure is isotropic 
at some given time to . Applying the change between Eu- 
lerian and Lagrangian coordinates = ^f-^~ twice 
- while neglecting higher orders - one obtains 



d 2 P ex n ox m d 2 P 



AkiAik 



dxidx; 



dxi dxj dXndXr, 



C, 



kk 



(•5) 



where C = exp(r^l) exp(r^4. T ) is the short-time Cauchy- 
Green tensor resulting from the integration of the flow 
map. The time r is chosen of the order of the Kolmogorov 
time so that it is possible to linearize the integration from 
diUj to diXj. 
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FIG. 1: The PDF of the gradient of the passive scalar V0 
using the Recent Fluid Deformation approach (top) and ob- 
tained from direct numerical simulations (bottom). 



FIG. 2: The PDF of the magnetic field B using the Recent 
Fluid Deformation approach (top) and obtained from direct 
numerical simulations (bottom). 



Applying a similar reasoning to the viscous term re- 

sults in the approximation i/AAij — — -^-Aij, T be- 
ing a characteristic friction time such that the Reynolds 
number 1Z is proportional to (^) [l|. In all simulations 
described below, the parameters T — 1 and r = 0.05 
were chosen which result in a Taylor-Reynolds number 
of 72. = 77. This value is in the range of the performed 
DNS simulations (see Table HJ. 

Finally, one obtains a stochastic ordinary differential 
equation (SDE) 



dAij — ( AikAkj —1 



^kk a 

3T A 



dt + dW, 



ij 



(6) 

where the stochastic forcing term W is added to repre- 
sent large scale forcing. The SDE ([6]) can be integrated 
with standard methods for stochastic equations UH)]. 



III. EXTENSION TO PASSIVE ADMIXTURES 

In this section we apply the Recent Fluid Deformation 
theory to the case of an admixture 9 which is passively 
advected by the fluid according to 



d t 9 + u • V0 = kA6 



(7) 



In order to study the dynamics of the gradient of the pas- 
sive field 8 we take the derivative of eqn. ([7]) and obtain 

dtidiO) + d t (u k d k 9) = Kdkkid.,9) . (8) 

Using again the Lagrangian material derivative ^ = 
(d t + Ujdj) eqn. © reduces to 

d 



dt 



(V0)i 



A ik (ve) k + Kd kk (V6)i 



(9) 



In this equation, the gradient tensor Aij is obtained from 
integrating eqn. ©. Modeling the viscous term as for 
eqn. and adding a stochastic forcing V, we obtain 
the stochastic differential equation (SDE) 



d(V0)< 



— (V0)i 



dt + dVi (10) 



This equation is integrated together with the equation ([6]) 
for the velocity gradient tensor. The forcing is chosen to 
be vectorial gaussian noise, delta-correlated in time , as 
for the velocity gradient. The forcing amplitude is V2di 
as for a standard random walk. The parameters for the 
result shown here are T@ = 1 and r = 0.05. 

In order to test the applicability of the Recent Fluid 
Deformation approximation, we performed direct numer- 
ical simulations for passive scalar turbulence. For this, we 
utilized the pseudo-spectral simulation framework LaTu 
9]. Table U summarizes the relevant parameters of the 
simulation. 

This model extension has been analyzed in [Io| in great 
detail, considering important phenomena like amplifica- 
tion of the gradient's norm and alignment to strain prin- 
cipal axes. But in contrast to this publication, we focus 
on the comparison of different turbulent systems. 

The top of Figure Q] shows the probability distribution 
function (PDF) of fluctuations of V6* obtained from inte- 
grating cqns. ([5]) and pU|) . In the bottom of Figure [1] the 
PDF obtained from the direct numerical spectral simu- 
lation is shown. The qualitative agreement is excellent 
and the tails of the PDFs of VQ approximately follow 
an exponential decay behavior as predicted by Shraiman 
and Siggia [11| for a Gaussian velocity field. 
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N 3 




»rm. 






V 


dx 


V 


L 


passive scalar 


256 3 


106 


0.0321 


0.0015 


4 ) 


< 10~ 4 


0.0245 


1.43 x 10~ 2 


2.02 


MHD 


256 3 


68 


0.41 


0.047 


2 ) 


< 10~ 3 


0.0245 


0.02 


1.47 



TABLE I: Parameters of the numerical simulations. N s : number of collocation points, = \/l5u rms -t//^: Taylor- Reynolds 
number, u rms : root-mean-square velocity, ek: mean kinetic energy dissipation rate, v. kinematic viscosity, dx: grid-spacing, 
n — (^ 3 /ek) 1/ ^ 4 : Kolmogorov dissipation length scale, L = (2/3E) 3 ^ 2 /e^: integral scale. Schmidt number and magnetic Prandtl 
number are equal to 1. 



IV. EXTENSION TO THE MHD EQUATIONS 

We start with the incompressible MHD equations in 
the kinematic regime neglecting the back-reaction of 
magnetic fluctuations on the velocity field. The equa- 
tion for the time evolution of the magnetic field B reads 



<9 t B + u • VB = B • Vu + 77AB 



(11) 



where the velocity field obeys the Navier- Stokes equa- 
tions In the Lagrangian frame this equation takes 
the form 



—Bi = B k A ki + rjAB . 
at 



(12) 



Again, using the same Recent Fluid Deformation approx- 
imation for the resistive term and adding a stochastic 
forcing, which is as above gaussian and white-in-time, 
we obtain 



dBi 



—BkAki 



32V 



Bi)dt + dUi 



(13) 



Note that this equation differs from (ITOl) by the trans- 
pose in A and - more important - the sign of the first term 
on the right hand side. Thus directions, where the mag- 
netic field is stretched, are directions, where the gradient 
of the passive scalar is compressed and vice versa. This 
very different influence of the stretching term results in 
a qualitatively different PDF for the magnetic field fluc- 
tuations compared to the V# fluctuations. In order to 
test whether the Recent Fluid Deformation approach is 
able to reproduce the different role of the stretching term, 
we performed direct numerical simulations of kinematic 
MHD turbulence. Here, we again performed simulations 
using the framework LaTu and the relevant parameters 
are given in Table HI The simulation was started with 
random initial condition for the magnetic field B and 
an already fully turbulent velocity field. External forc- 
ing was only applied to the velocity field by keeping the 
large scale Fourier-modes k < 2 constant. The top of 
Figure [2] shows the PDF of magnetic field fluctuations 
calculated using equations (J5]) and (fT3"|). This PDF has 
to be compared to the PDF obtained from direct numer- 
ical simulation shown in the bottom of Figure [2] It is 
again remarkable that the simple stochastic ODE model 
(fTB")) is able to capture correctly the shape of the PDF of 
magnetic field fluctuations. 



Since we considered the kinematic MHD equations in 
the Dynamo regime, the PDF of the stochastic model 
could not be obtained from a time series as in the case 
of the passive scalar since the magnetic field is growing 
exponentially. 

To obtain meaningful statistics ensemble averaging was 
used instead of time sampling. Numerically, a random 
initial state was generated and integrated a sufficiently 
large number of time steps, to retain no statistical in- 
fluence of the initial state. The final state enters the 
statistics and a new realization is generated via another 
random initial condition. The PDF was obtained by sam- 
pling over 10 6 initial conditions. 

For illustration of the kinematic dynamo effect and a 
consistency check with direct numerical simulations of 
the kinematic MHD equations we tracked the local mag- 
netic energy B 2 (X(i),t) in Figure [3] for Tg = 6, a value 
well in the dynamo regime of the modell. For this sim- 
ulation only the velocity gradient itself was driven by 
stochastic forcing, as outlined above, while the magnetic 
field was undriven. The forcing amplitude was scaled 
with V2d£, analog to the scaling for a standard random 
walk. 

The parameter r of the model was chosen to be r = 0.05. 
The magnetic diffusion time Tg scale is varied to explore 
the reaction of the dynamo effect. Below the critical pa- 
rameter of approximately Tg ~ 4 we observed no growing 
magnetic energy. 

The corresponding growth rate, estimated from an ex- 
ponential growth ss exp(7£/r), is about 7 ~ 0.032. The 
growth rate compares well with the growth rate of our 
DNS simulation, where we obtain 7 ~ 0.03. Although 
the identification of the correct Kolmogorov timescale r 
is not without uncertainties, it is remarkable that this 
value of the growth rate is in the range of data obtained 
by kinematic dynamo simulations (see also Table I in 12 1 
and references therein). 



V. CONCLUSION & OUTLOOK 

In this paper we have shown that the natural extension 
of the Recent Fluid Deformation model for Navier-Stokcs 
turbulence [l[ to the case of fluctuations of the gradient 
of a passive scalar and to magnetic field fluctuations is 
able to produce probability distribution functions that 
agree well with PDFs obtained from direct numerical 
simulations. The PDFs from the stochastic ODEs ([51 



4 



[TOl H~3l) are obtained with a fraction of the computing re- 
sources necessary for direct numerical simulations. The 
next step is to include the back-reaction of the magnetic 
field on the fluid flow. If this can be managed, then gen- 
erating magnetic field fluctuations using the Recent Fluid 
Deformation model is a tempting alternative for e.g. 
the problem of cosmic ray propagation where magnetic 
field fluctuations are generated by other means which 
are not able to capture intermittency effects (see [HI, [l4[ ) ■ 
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FIG. 3: (Color online) Temporal growth of the magnetic field 
energy. To be well in the dynamo regime we chose Tb = 6. 
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